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ABSTRACT 

The Mavstar Global Positioning System Is a highly accurate, space based navigation systen providing all veather, 24 
hour-a-day service to both military and civilian users. The lavstar systea provides a Gaussian position solution with 
four satellites, each providing Its epheMrls and clock offset vith respect to GPS tine. Currently, GPS Is building 
towards the full 24 satellite constellation. 

The GPS Master Control Station (MCS) Is charged with tracking each lavstar spacecraft and precisely defining the 
ephenerls and clock parameters for upload Into the vehicle's navigation nessage. This paper briefly describes the 
lavstar systen and the Kalnan Filter estimation process used by the MCS to determine, predict, and quality control each 
satellite's ephenerls and clock states. Soutine performance is shown. Kalnan Filter reaction and response Is discussed 
for anomalous clock behavior and trajectory perturbations. 

Particular attention Is given to MCS efforts to Improve orbital adjust nodellng. The satellite out of service tine 
due to orbital naneuverlng has been reduced in the past year from four days to under twelve hours. The planning, 
reference trajectory nodel, and Kalnan Filter aanageient improvements are explained. 

Finally, this paper will summarize the future work to iaprove Kalnan Filter and orbital adjust performance. 
Keccomendations will be made for future systems requiring precise ephenerldes. 


IMTRODDCTIOM 

The Mavstar Global Positioning Systen (GPS) Is a satellite based radlonavlgatlon systen providing worldwide all 
weather coverage to both civilian and military users. The Mavstar systen promises a revolution in all activities 
requiring precise navigation or positioning. GPS signals will provide the precise positioning service to authorized 
userB of 16 meters spherical error probable (SEP) and IN seter circular error probable (CEP) to standard positioning 
service customers. The standard positioning service Is subject to change according to Baited States national 
interests. GPS is made up of three segments. The space segment consists of the orbiting satellites and provides L-Band 
signals with nodulated data to the world. The User segnent represent* the customers who receive and utilize the 
navigation data. Finally, the control segnent comprises a systea of L-Band monitor stations, S-band ground antennas, and 
a control center to monitor the satellites health and periodically upload new navigation parameters. GPS is currently 
the largest satellite constellation dedicated to a single purpose, and It's still growing to the planned 24 satellite 
constellation consisting of 21 satellites and three on-orbit spares.. The Master Control Station (MCS), located at Falcon 
Air Force Base Is responsible for the nalntenance of the GPS satellites and their payload. The navigation payload of 
each satellite is regularly updated with orbital paraneters, atomic frequency standard states, and aiiwi c data for 
broadcast to the user community. This paper will briefly describe the Mavstar system, the arrangement and organisation 
of the MCS, and the Kalnan Filter algorithms used by the MCS. The paper will go on to describe the ways In which the MCS 
manages those algorithms. Nominal performance as well as performance under anomalous conltlons will be discussed. 

The improvement in orbital adjust modeling will be explained. 


THE XAVSTAK SYSTEM 

The GPS satellites operate at an altitude of approximately 21183 Km In a near circular orbit with a 12 hour 
semisynchronous orbit. Orbit inclination for the operational satellites is 5J degrees. Each satellite transmits two 
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L-Baod signals. The LI signal Is a radio carrier frequency of 1575.42 MH*. The LI signal Is nodulated by two psendo 
randan noise (P1H) codes known as the course acquisition (C/A) and precision (P) codes which are in phase quadrature with 
each other. The P code is a PM strean of 11.23 Mbits per second conpared to the C/A code's rate of 11.23 Kbits per 
second. The 12 carrier operates at a frequency of 1227.6 MHz and usually only contains the P code. Both the P and C/A 
PIH codes f p nt al" the nodulated data representing the navigation nessage fron the transacting satellite. This 
navigation nessage is transferred at a rate of 5» bits per second. The navigation nessage contains the satellite’s clock 
offsets fron GPS tine, the satellite's precise epheneris, and course alaanac lnfornatlon for the entire operational GPS 
constellation. Other paraneters are included in the navigation nessage to indicate the health of the GPS satellite, give 
single frequency users ionospheric correction lnfornatlon, and provide Universal Tine, flsers can either receive the 
single C/A code, or use the C/A code to acquire the P codes on LI and L2. All clocks in the Mavstar Systen are 


synchronized to GPS tine. 

Once the user reads the navigation nessage fron each of the visible tracking satellites, the user can synchronize 
his tine to GPS tine, calculate the apparent range to the satellite, and solve for position using the fornula: 


Where, 


PR! 


- V (X i • X u ) 

<x i ,Y l ,z l )T * 

(X , Y ,Z ) T • 
u u u 


2 * (Vi - V 2 ♦ Ui - z u )2 4 C ^U 

Satellite Inertial Position 
User Inertial Position 

lonospherically Corrected User Pseudorange Measurement 
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( 1 ) 


c - The Speed of Light 

- User Clock Offset from GPS Time 

Since the satellite inertial positions at a given GPS tine are derived fron the navigation nessage and pseudorange 
neasured fron the PM code tine differences, the only unknowns are the user’s inertial position and the user's clock 
actual offset fron GPS tine. Since their are four unknowns, four pseudorange neasurenents fron four different satellites 
are required to find the user's inertial position and clock offset. Table I shows the specifications for GPS user 
accuracy conpared with other navigation systens. Obviously for this systen to work properly, knowledge of the 
satellite's epheieris and clock offsets with respect to GPS tine nust be precise. The GPS Operational Control Segnent 
(OCS) is charged with keeping the navigation nessage accurate and the spacecraft in good operating condition. 

The Operational control segnent consists of five L-Band Monitor Stations (MS), three S-Band Ground Antennas (with 
two additional antennas available on a part tine basis), and the Master Control Station (HCS). The locations and 
coverage of these ground stations is shown in Figure I. The L*Band tracking data is transmitted fron the monitor 
stations to the HCS. The MCS uses this information to generate new orbital elements and clock states which are 
transmitted through the ground antennas to the satellite. Information on the satellite's status is also gained fron the 
ground S-Band telemetry. This system is pictured in Figure II. The GPS MCS is the hub of all the activity in 
the coaand and control of the Mavstar systen. 


THE MASTER CONTROL STATION MANNING 

The GPS Master Control Station's operations center is staffed with seven members of Air Force Space Comand's Second 
Satellllte Control Squadron. This group of seven people comprise the basic space operations crew. There are five crews, 
each working for a period of six days followed by four idle days. The off days provide for crew rest, as well as training 
and standby days. These operations crews maintain the Mavstar constellation's state of health and navigation performance 

around the clock, every day of the year. ...... 

There are currently six positions in the operations center filled by the seven crew members. The Ground Controller 
is responsible for the comranicatlon links between the MCS and each of the four ground antennas and five monitor 
stations. The Satellite Operations Officer is the position tasked with making contact with a given satellite from the 
MCS through a ground antenna to monitor a given vehicle's state of health and transmit any necessary commands. The 
Satellite Engineering Officer is responsible for ensuring that a spacecraft is in good operating condition, all commands 
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TABLE I: 


Comparison Of GPS Specifications 
Positioning Service. 


To Other Navigation Systems. 


Accuracies Are Given For The Precise 
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Position accuracy lc degraded mainly heron* 
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90 minute*. For in* in How moving vehicles 
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dual frequency measurements. 
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2 SNU-S4-1 Sgcifiotion for USAF Standard Form Fit and Function 

•Mium Accuracy Inertial Navigation SeSunit. AcufaarlUi — ' 

X Ot pendant on integration concept and platform dynamic* 


GA VISIBILITIES (5* ELEVATION) 


MS VISIBILITIES 




West LONGITUDE Ea«t 


COSPM - Colorado Springs Monitor Station 
ASCNM - Ascension Island Monitor Station 
DIEGOM - Diego Garcia Monitor Station 
KVA TM — Kvajnleln Atoll Monitor Station 
HAWAIM - Hawaii Monitor Station 

States on a part-time basis. 


West 


LONCITUDE 


East 


CTS - Colorado Springs Tracking Station 
CAPEG - Cape Canaveral Cround Antenna 
ASCNC - Ascension Island Ground Antenna 
DTEGOG • Diego Garcia Ground Antenna 
KWAJC - Kvajalein Atoll Ground Antenna 


NOTE: CTS and CAPEC give coverage over the Continental United 


FIGURE I: The CPS Ground Antenna and Monitor Station Locations and Their Respective Visibilities 
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sent to the spacecraft are valid, and to respond to any anomalous behavior. J 

The Satellite Analyiis Officer wmitors the l-Band tracking data fro« each of the aonitor stations and verifies the 
navigation performance of the GPS constellation. The Deputy Colander takes care of all the bookkeeping and upchannel 
reporting of the crew’s activities required in an Air Force operation. Finally, there's the Crew Coaander who sanages 
the entire crew and perfons real-tlie rescheduling of the operations center's activities. 

At the present time, there are two Satellite Operations Officers while every other position is filled by one crew 
seiber. The Ground Controller is at this tine the only enlisted person on crew. Every other position is staffed by 
junior officers. The next ten years will probably see the Satellite Operations Officer positions grow to three lots and 
transition to enlisted personnel. The Deputy Cosnander position is also likely to be taken over by an enlisted person. 

In addition to the operations crews, the Second Satellite Control Squadron supports these crews with several other 
sections. The Space Operations section directly supports the operations crews. The Scheduling Section produces a daily 
schedule of events which the crew follows. The Training and Standardiiatlon/Evaluation sections guarantee crew 
proficiency in their Jobs. There’s also an engineering support group which is broken down into sections responsible for 
the analysis of the spacecraft bus. navigation payload, ground system, and systen database. Finally therein a section 
dedicated to interfacing with various users and a Comand Section. All in all. there are approximately IN people in the 
Second Satellite Control Squadron responsible for the operation of the GPS constellation. 

It's the Satellite Analysis Officer crew position and the navigation payload section of the engineering analysis 
group which work with the contents of this paper. This group of people manages the Kalman Filter algorithms and results, 
guaranteeing the end product to the GPS positioning/tining users is usable, accurate, and within all published 
specifications. 


GPS OPERATIONAL CONTROL SYSTEM 



c=> NAVIGATION SUPPORT 

CONSTELLATION MAINTENANCE 


FIGURE XI: The CPS Operational Control Segment. L-Band signals transmitted by each satellite are tracked 

by the Monitor Stations. Navigation Uploads, Telemetry, and Commands are routed through the 
Ground Antennas. The Master Control Station is the central processing location. 
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THE MGS KALMAN FILTES P10CESS 


The tiro fundamental measurements made by the Mavstar systei are pseudorange (PR) and accumulated delta range (ADI). 
Pseudorange measures the distance between a satellite and a monitor station on the Earth by accounting for the total time 
between transmission and reception. This measurement is a "pseudo’ range rather than a "true" range as a result of 
delays encountered on the signal's path to the Earth's surface. The equation used for pseudorange by the MCS is: 


Where, 




Pseudorange Measurement 
Speed of Light 

Time Signal Received by Monitor Station 
Time Signal Transmitted By Satellite 


( 2 ) 


The transmission time is represented by an integer number of r-counts where one i-count is equal to 1.5 seconds. The 
satellite aid toil tor station tines are referenced to GPS tine for synchronization. 

,.v. J elt f ran9 f JJJ 1 ] } 8 1 ““wement of total phase accumulated between the transmlslon and reception of 

the carrier signal. Accuwlated delta range Is gl?en by the expression: ^ 


ADR 
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( 3 ) 


Where, 


ADR 

f 



Accumulated Delta Range 

Carrier Frequency 

Monitor Station Phase Measurement 

Transmitted Satellite Phase 


The only use for ADR's within the GPS control sepent Is to aid in the smoothing of pseudorange data. ADR's are 
available primarily for the GPS user coumunlty to calculate relative velocities. 

. J* *“. ter Con f ro1 s ^ tloa continuously tracks and collects pseeudorange data from each satellite visible to the 
■onitor stations. A pseudorange measurement Is made every r-count on both the 11 and 12 signals received Data 

at llMte lntem1 *- ^ ton in Figure III, tbe MCS begins to store data at each fifteen 
■lnute point. Over each fifteen minute interval the satellite transmissions art checked for correct parity and C/A to P 
code Handover Word. Each measurement nst also be within limits for slgnal/nolse ratio, code and phase slips, code and 
carrier lock, and first and second difference tests. At tbe end of tbe fifteen minute interval there is a five minute 
waiting period. This waiting period Is nsed to retrlewe any data which may reside In a monitor station's data buffer due 
to short communications outages. After the five minute wait, the data which passed all editing crltertla for raw 
pseudorange measurements daring the fifteen minute Interval are now ready for processing. 

Each LI and L2 measurement passing the editing criteria is now used to detenlne the lonospherlcally corrected 
pseudorange. The layers of the Earth's atmosphere approximately between (« and (41 Km altitude constitute tbe 
ionosphere. Tbe ionosphere can vary greatly over the surface of the Earth for such reasons as solar activity, the 
effects of man, and the variation between day and night. The ionosphere affects GPS by bending the transient carrier 
signals. Fortunately, the L-Band signal frequencies nsed bt the Mavstar system exhibit the property that the Index of 
refraction is proportional to the Inverse square of the frequency. Since the two different carrier freqnencles, LI and 
L2, measure the same distance between the satellite and the monitor station, the proportionality constant can be 
determined independent of the actnal ionospheric conditions or the elevation of the satellites. The proporlonalltv 
constant Is given by: 
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where 



PR 1 * Direct LX pseudorange measurement 
PR 2 * Direct L 2 pseudorange measurement 
* LI carrier frequency • 1575*^2 MHz 
f 2 * L2 carrier frequency * 1227.60 MHz 

Using the LI neanreaent, the ionospheric tine delay nov becones: 



k t 


2cf 


2 

1 


( 4 ) 


(5) 


So each pair of LI and L2 pseudorange neasnrenents received at every z-connt is used to transform the raw pseudoranges to 
an ionospherlcally corrected psendorange with the equation: 


p “i 



( 6 ) 


The ionospherlcally corrected psendorange is the only calculation which occurs to the signal before filter 
processing. Other phenoaenon which affect the signal, such as troposherlc delay, relativity, and the free space delay 
are accounted for within the filtering systen. 

The ionospherlcally corrected pseudoranges are now snoothed over the IS ninute collection interval to reduce the 
neasurenent noise before Kalian Filter processing. A polynomial fit is node to the ionospherlcally corrected 
pseudoranges for each satelllte*nonltor station neasurenent pair. The snoothed pseudorange is detenlned directly from 
the polynomial and is set at the beginning of the fifteen ninute Interval. 

The vector containing snoothed pseudoranges for each satelllte-aonltor station pair is defined as the observation 
vector. This observation vector will be processed by the MCS Kalian Filter to produce a state vector for each satellite 
and aonltor station as shown in Figure IV. The state vector for each satellite Includes its inertial position and 
velocity, a scaling parameter (Kl) and an acceleration parameter (K2) for a solar pressure force model, and clock bias, 
drift, and drift rate. The aonltor station state vector Includes clock bias and drift, and tropospheric height. 

To reduce the aagnltude of the estinatlon problem, the MCS uses an Independent partitioning scbeae. Each partition 
is aade of state vectors for up to six spacecraft and Includes all monitor station states. This scheme reduces the 
processing load of the MCS, allows for isolation of satellites with poor performance, and gives the operator insight into 
the ground stations' states between partitions. 


00 



A - 15 Minute Data Collection 

Period (600 PR Measurements) 

B - 5 Minute Wait Period 

C - Approximately 2 Minute 
Processing Interval 


FIGURE III: 


Clock Representing The Time Intc'rvals 
Data. The Cycle Shown Above Starts Every 15 Minutes. 


Ir Which The MCS Collects And Processes Psemlor«n*e 
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12 


SPR 
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SPR 


ij 


Jf • Observation Vector 
SPR - Smoothed Paaud orange Kt«»ur«m*nt 
1 • Satellite Number 
j • Hon i tor Station Rubier 


MASTER CONTROL STATION 
KALMAN FILTER 
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FICl'RE IV: 


The MCS Ka loan Filter Uses Smoothed Pseudoranges To Determine The Satellite And Monitor Station States 
Shown Above . 


The purpose of the MCS Xalaan Filter Is to estimate the corrent states In a partition given the available 
aeasurenents in the observation vector. In the MCS aanlfestatlon of the Xalaan Filter, the current state process is 

u l r ? 6 ?! Ct J 0 / refereace 8?8t “ “* then related to an epoch state residual process. The observation 
process is llnearited about an apriori estlaate of the state vector and likewise related to to the epoch state residual 
process. The foraulation of the MCS Xalaan Filter is described in the following paragraphs to give soae context to the 
subsequent discussion of GPS Kalian Filter operations. 

. f In J h ®. lB Pfeaentation of the Xalaan Filter, a predicted reference trajectory and a set of reference clock states 
exist. A staple relation is drawn between each state and the reference states at a specified tine. For exaaple the 
current state residual is defined as the difference between the current state and the reference state at a current tiae: 


iref(t) 


( 7 ) 


epocb T tiar° Ch ***** re8ld “ l is defined 88 the dlfferen « ktveen the epoch state and the reference epoch state at the 



&o ~ &ref(o) 


( 8 ) 


The epoch states are napped to the corrent states through a nonlinear state propagation aodel (F| . . . 
reference state at current tiae can be f ound froa the reference epoch state at epoch tine: 


I) such that the 


*ref " F,t * V *ref (0) ^ 


( 9 ) 


The current state can be found froa the epoch state and their associated tines with: 


it 


Ft*. 
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The relationship between the state residuals and their respective napping functions is pictured in Figure V. Mote 
that the reference state process shown in Equation 9 Is a completely deterainlstlc function. 

For every state estinating partition, the state process is noddled in discrete tine and given as: 

x(T k+l ) - F|i k+l . v X(T k )| + w(i k ) (ll) 

Where, Ji - State Vector 

v ■ State Process Noise 

■ Time of Current Data Interval 

T k+1 « Time of Next Data Interval 


The state and state process noise vectors are dimensioned by the nunber of states in a partition. This current 
state process is linearized by substituting Equation 7. at the discrete tine, into Equation 11 and then performing a 
first degree Taylor Series expansion on the current state residual ten. With the substitution of Equation a 
discrete tine into the expansion, the current state residual is found to be: 




dF(T k+1 , y X 1 


dX 


*<V + «(x k ) 


( 12 ) 


X-x 


ref 



Mapping Functions. 
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By following tbe me logic, the current state residual is related to the epoch 
napping function with respect to the reference epoch state. Substituting Equation 8 
expansion on the epoch state residual tern, and substituting Equation 1 results with 


>x( V 


FfT. . T-. X 


ax 
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state residual by linearizing the 
into II, perforning a Taylor Series 
the equation: 


(13) 


X-2 0 


If the two discrete tine forns of Equation 13 are substituted back into Equation 12 and algebraically nanlpulated the 
final result is the epoch state residual process used by the MCS: 




WV 




-1 

s<V 


(14) 


The state process noise is assuned to be a white gausslan process with zero nean such that the expected values of 
the state process noise have the font: 


Elw(T k )J - 0 and ] - Q(X k ) ( , ; 

5^ • Kronecker Delta 

« Any Other Data Interval 

where Q represents the covariance Matrix for each state. The covariances have a constant value defined in tbe MCS 
database, and are added at each Kalnan Filter data Interval. 


The observation process is aodelled in discrete tine with the nonlinear equation: 




h( 5 t J 


(16) 


The observation vector and process noise are dlnensioned by the total nunber of observations nade by a partition. 
This equation relates the current state vector to the snoothed pseudoranges with the state transition Matrix, h( 1 with 
the addition of Measurenent process noise. This observation process is linearized about an apriorl state estlnate and 
then expressed in terns of the epoch state residual. The resulting linearized equation is: 
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As show in Equation 2 the pseudorange aeasures the tine delay between the tine a nonitor station received the given 
signal. This tine difference can be broken down further to its conponents such that the observation nodel takes the 

torn: 


h(X(T k )) - cj 

’ AvV 
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Where* 

A T t ’ 


A t fs " 
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(18) 


FS 

Tropospheric Delay 
Free Space Delay 

Relativistic Delay 
Monitor Station Clock Offset 
Satellite Clock Offset 

Signal Transmit Time * Data Interval Time 
T r - Signal Reception Time 

The last three terns inside the brackets give the tiling offsets between the nonitor station and satellite clocks. The 
first two terns represent the propogatlon delay between the satellite and nonitor station. Renenber that the ionospheric 
delay has already been accounted for in the dual frequency corrections nade before snoothlng. The clock offsets are 
slnply the tine bias between either the satellite or nonitor station clock and GPS tine. The free space, tropospheric 
and relativistic delays will briefly be described. 

The free space delay Is experienced by a satellite's transnltted signal due to the physical propogatlon path length 
in vacuun. The free space delay is given as: 

r ms (t r*I 


Awv V 


*sv<V 


(19) 


a Satellite Inertial Range 


The nonitor station reception tine and the position vector for the nonitor station phase center are found through an 
iterative process given the known transit tine. 

The troposphere is the portion of the atnosphere which fills the volune between the Earth's surface to an altitude of 
about 12 kiloneters. The najor effect of the troposphere is to give an apparent increase in the signal path length. To 
help in the evaluation of tropospheric delay, three atnospheric paraneters are neasured at the nonitor station recieving 
the signal. These paraneters are atnospheric pressure, tenperature, and dew point. Since different locations around the 
world can have widely varying conditions of the troposphere, the tropospheric height is nodeled as a state variable in 
the estinatlon process of the Kalnan filter. The tropospheric height is the actual altitude of the troposphere at the 
nonitor station. The tropospheric delay is found with the equation: 
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Monitor station Barometric pressure converted to 
K i lopasca I s 

Monitor station measured temperature converted to "the 
Kelvin scale 

The r-adiai distance from the Earth's center to the 
meteorological sensors. 

The radial distance from the Earth's center to the 
monitor station antenna 
Troposohertc dry radius 
Sate M 1 te elevation 
Tropospheric neiynt 
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= Assumed integration constant 
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ftB ru SS Ut J V rMt J!l 11 k ! XperleBC ! th * effMt 0f ,pecU1 relatl?lt ? dfle t0 it* speed with respect to an observer 
on the Enrth. To the ground observer, the satellite clock will appear to be behind. The sane satellite <nii «ito f,Ii 

nlf f J^ni 1 rehtiTlt L be “ Me of the orblt ’ s lwer 9»»itational potential. To the sane ground observer the 
te l V ?? ¥l11 app< ff ahead ‘ ** effect of 9 *wral relativity Is greater so to compensate, the satellite's 

inhibit TJ 1? mher th “ the 11,23 10,1 of tbe >0Blt0r 8tatlon docks. Wth the satellite 

In orbit the two clocks will appear to be synchronized* 

,, cloc * ^hind the ,0ttltor *t*tlon clock would completely account for the relativistic effects 

if i m “ ? lllte * orbit ! ere P* rfectl T circular. Since the satellite Is never In a perfectly circular orbit a periodic 
relativistic variance needs to be included with the equation- 1 periodic 




sin(E) 


( 21 ) 


whert 

e * Oroit eccentricity 

fj =* Universal gravitational paranefer. 
a - Sem i-(T3jor axis Of orbit 
c 3 The ipee j of ' i ght 

E = Eccuni r ic Anomaly From Keple r * a Equal ion 


th* iwlif,!?, l eac i dat \ co ! 1 ! ctlon lnter * 41 ' 4 Kall “ Filter solution Is used to find the current state vector with 
fw 1 ??!? !j U fj ®P° cb 6t4tc P ncest (Equation 14) «ad obsemtlon process (Equation 17). A tine update occurs 1 
first. In the tine update, the aposterlorl state estlnates fron the previous data interval hecone the aprlori state 
estlnates for the current data interval. Process noise is added to the aposterlorl covariance natrix to define the 
aprlori covariance natrix for the current data interval. The equations dictating the tine update are: 


■ Saw 

( 22 ) 

!o ( w - Vv * ( B<T k ) ] Q(T k ) [ B<T k ) ] T 

“ Apriorl estimate of epoch state residuals given past observations. 

■ Aposterlorl estimate of epoch state residuals given past and current 
observations. 

* Covariances of the epoch state residuals. 

* Process Noise 

■ [0(T k+ l ,T O ) ] [ G ^ T k^] wherc; [ G ^ T k{] ls a wie * htln S matrix. 

After the tiwe update occurs, a Kalian Gain is calculated wing the aprlori states and expected covariance latrlx: 

$ “ to [»<&> ] T [(wXtifio l H( 5t>] * M J" 1 <23 > 

The aposterlorl estiuate of the epoch state residuals for the current data interval is now found fron the equation: 

<5*o • 4[* t - h <5 t )] (24) 

to -Jo - * !«<&>]& 

Using Equations 8 and II, the current states can be found fron the epoch state residuals. 

With the linearized epoch state process, if the current state deviation fron the reference states becomes too 
large, the quality of the MCS Kalman Filter solution can degrade. To avoid this problem, the ACS updates reference 
trajectories every four weeks and whenever linearity Units are exceeded. 


Where, 


t>W 

$0<V 

.w 

Q(T k ) 

[ B( V] 
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DAYS SINCE REFERENCE TRAJECTORY EPOCH 


FIGURE VI-A: Typical Epheraeris State Performance. Inertial Position Offsets Have Been 

Into Body Radial. Along-Track. and Cross Track Deviations. The Offsets t 
Epoch are Zero. On the 28th Day. The Current Reference Trajectory is Replaced with a 
New Reference Trajectory Based on the Current State Estimates at the New Epocn. 
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FIGURE VI-B: 


Satellite Cesium Atomic Clock State Bias and Drift Offsets From GPS Time. These Plots 
are also Typical of Normal Performance. The Reference States for Clocks in the CPS 

"ffast!r*c? D w.* r “ lne r fro " G P S T1 "»- GPS Time Is Derived from Either a Single Designated 
Master Clock or a Composition of All the Clocks in the CPS System. Note That There 
lsn t any Apparent Influence on Clock States From Ephemerls States or Vice Versa. 


j MCS KALMAN FILTER OPERATIONS 

| Figure VI-A Indicates the typical performance of current state residuals over a period of one week, approximately 

three weeks after a reference epoch. Aging of the reference trajectory states is evident over this period. The aging 
j occurs because only najor orbital perturbations are modeled by the Master Control Station. Figure VI-B shows typical 
cesium atomic clock performance compared to GPS time over the same period as the ephemerls states. 

Epoch State Residuals, Current State Residuals, Epoch States, and Current States, as well as plots over a twenty- 
four hour period are all available for viewing to the Satellite Analysis Officer (SAO) at the end of each data interval 
Though always available to the SAO, this information is usually only accessed when additional data on a given situation 
is needed. 
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The difference between the ionospherically corrected pseudorange and the pseudorange calculated from the apriori 
state estimates Is called the pseudorange residual. The SAO observes the pseudorange residuals at the end of each data 
processing Interval. The value of the pseudorange residual is usually below three neters. Values inch above six meters 
are usually the first, and best indicator that there is a monitor station proble*, satellite problen, or possibly a poor 
Kalian Filter iodel. If the pseudorange residual is above a dynaiic liiit based on the covariance latrix, the 
neasureient obtained during the data interval is not used in the update of the Kalian Filter model. 

When there are probleis with the state estimates of a satellite or lonltor station, the SAO has rciatively few 
options to reiedy the given situation. One option is to allow the process noise to add over tine such that the 
neasureient falls within Units sole tiie in the future. This option generally isn't acceptab e due to elther the to a 
tine involved or the size of the rejected neasureient. If the error can be attributed to a satellite or wnitor station 
clock state the SAO is capable of specifically nodifying the clock state based on the observed pseudorange residuals. 
Another option left to the SAO is to change the value of the covariance latrix IP] in Equation U. This action has been 
inappropriately called the "Q-Bump". It’s possible to lodlfy either nonitor station, satellite, clock, epheieris, or 
lolS pressure T state covariances. Any subset of these state covariances can be "O-Baped" the “*e time. The value 
of the "Q-Buiped" covariances is prescribed in the MCS database and can vary with different situations. It s also 
possible to reinitialize the estination proble*, but this action is usually a drastic step. In soae situations, various 
MCS database parameters can be changed which affect the behavior of the Kalian Filter algorithms. 

There aretwo other tools used by the SAO. The first tool is the capability to usk the Kalian Filter. Masking 
prevents the update of the Kalian Filter state estliates with new data, thereby forcing only a tiie update to occur each 
data interval. The second tool is the capability, on a non-interference basis, to reprocess the Kalian Filter over a 
twenty four hour period with varying conditions. This reprocessing capability is quite powerful and chiefly responsible 

for the resolution of nany probleis. 

The navigation data providing the CPS user connunity the precise ephenerides needed for accurate positioning are 
based on the Kalian Filter State est Hates. When a navigation data upload is generated for .transmission to a Mavstar 
spacecraft, the satellite’s epheieris is predicted by using the aposteriorl state estimates to differentially correct 
the reference trajectory. The clock states are predicted aerely by propagating the aposteriorl clock state estimates 
into the future. These epheieris and clock state predictions are transmitted to the spacecraft with S-Band and, through 
the satellite, transiitted to the user community l-Band at the appropriate tiie. ...... 

The Master Control Station's insight into the quality of the upload prediction is performed with two statistic 1 
measures using the Kalian Filter state estimates. The first measure is the Estimated Range Deviation ERD). The ERD is 
the difference between the aposteriorl estimate of pseudorange and the pseudorange derived fro* the upload prediction 
for a particular vehicle. The ERD is calculated for a mathematically determined set of locations around the world and 
each GPS satellite visible to. those locations. Normally ERD’s grow over time indicating a slow degradation of the 
upload prediction. Figure VII shows this aging process. The need to upload each Navstar spacecraft with new navigation 
data on a dally basis is apparent in Figure VII to maintain the specified navigation performance. OD s usually vary a 
little more than shown in Figure VII due to such aspects as *onitor station visibility, timing of the upload creation, 
and clock movement. The second leasure of performance is the Observed Range Deviation (ORD). The 1 3RD is the difference 
the ionospherically corrected sioothed pseudorange and the aposteriorl pseudorange determined by the navigation upload 
state prediction. An ORD is produced for each satellite visible to a nonitor station. The ORD can be a good indicator 
of navigation performance. Should the ERD's and ORD's rise above a tolerance, set by the MCS as II Meters, a 
contingency navigation data upload will be accomplished to maintain user specified accuracy. Equations describing the 
Pseudorange residual, ERD, and ORD are given below: 


PRR 

- SPR tJ - 

PR(X) J . * Pseudorange Residual 

- ij 

ERD 

- 

*** ^upload predict ion^ij 

- Esc lasted Range Deviation 

ORD 

’ SPR ij - 

^^*upload predict ion^ ij 

Observed Range Deviation 


I 


(25) 
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FICURE VII: 


Typical Upload Prediction Performance a* Reported by ERD's. Performance la Measured 
gainst Range Based on Current State Estimates. Navigation Uploads Contain Eohemarls 
and Clock State Prediction. Baaed on the Most Current State Estimates! ? 


RESPONSE TO ANOMALOUS EVENTS 

In January 1989, a Cesium atoalc clock aboard Navstar vehicle number 11 failed. This failure necessitated a clock 
witch on board the spacecraft to a Rnbidlna atoalc clock. A "Q-Bump" on the satellite's clock states vas performed 
after the hardware failed. The Xalaan Filter estlaates had stabilised to new clock states In approximately six hours. 
This performance is shown in Figure VIII. The transient during this twelve hour period Is readily observed. One aonth 
prior to the clock witch, the spacecraft bus section of the 2SCS engineering support group, predicted the cesiua clock 
failure. As the tlae of the predicted failure approached the bad clock's performance was deteriorating. This 
deterioration is barely discernible In the clock drift plot shown in Figure VIII. Even with the deteriorated clock, 
performance of the satellite's navigation payload was maintained through a series of Kalman Filter clock state 
modifications and contingency navigation data uploads. On the fifth day of the year, maintaining performance through 
this means became impractical, so the satellite was set "Unhealthy" to the user community and the atonic clock was 
switched. At this time, the Kalman Filter was masked to prevent update of the satellite’s clock and ephemeris states 
and the new Rubidium clock was allowed to "warm up". On the sixth day of the year, the clock state "0-Bump" was 
performed. Clock performance was observed over the next several day6 and Navstar 11 was set back "Healthy", on the 
tenth day of the month. 

Figure IX illustrates the impact of a satellite trajectory perturbance. In this case, Navstar 18 came out of a 
solar eclipse with an attitude error. This error was sensed by the spacecraft's Attitude and Velocity Control Subsystem 
(AVCS), and compensated with an attitude thruster firing. The Impact on the trajectory In this instance was readily 
followed by the Kalman Filter State estimation process. In the situation shown In Figure II, only a contingency upload 
was required to naintain navigation performance and a new reference trajectory was built approximately one week after 
the event to prevent linearity failures. 

Satellite Vehicle thruster momentum dumps look very similar to Figure II, but may be more severe. Nonally the 
satellite will reduce reaction wheel momentum through variable setting magnets which torque the vehicle. Occasionally, 
momentum will saturate the satellite's reaction wheels at which time the vehicle's AVCS will command a thruster momentum 
dump. Sometimes this event may not be seen by the NCS due to visibility constraints. The momentum dump will manifest 
itself with growing ERD's and ORD's with respect to the ephemeris estimate and after twenty four hours, the pseudorange 
residuals will be rejected since the resulting states are above the calculated limits, indicating a mlsnodelllng 
problem. In this situation, the Kalman Filter Is reprocessed with an ephemeris covariance "0-Bump" near the expected 
time of the thruster momentum dump. Usually only a contingency upload 16 needed to maintain navigation performance 
followed up later with a new reference trajectory based on the changed state estlaates. If the pseudorange residuals 
ERD's and ORD's are over tolerable limits the satellite is set "Unhealthy" to the user coHunlty. 
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SATELLITE ALONG-TRACK 
DEVIATION 



FIGl’R£ IX: The Affect on Ephemerls State Estimates After an Attitude Thruster Firing on NAVSTAR 10. 

Satellite Thruster Momentum Dumps and Well-Modeled Delta-V Manuevers Show the Same 
Characteristics as this Plot* Radial and Cross-Track Deviations are not Shown for Brevity. 
There wasn't any Noticeable Clock State Perturbations due to this Ephemerls Change. 


DETA-V PLANNING AID MODEUG IMPROVEMENTS 

Prior to 1989, satellites undergoing Delta-V orbital maneuvers were usually 6et "Unhealthy* to the user coMunity 
for a period of four days after the trajectory change. The engineers and the analysts in the 2 SCS didn t believe that 
this type of performance was acceptable. There were several efforts made to shorten the total Unhealthy tine due to a 
Delta-V. 

The first improvements were Bade in the area of Delta-V planning. The location of a GPS satellite is required to 
be at a specific longitude of ascending node with a tolerance of 47- 2 degrees. When a satellite originally approached 
this tolerance boundary, a station nalntenance Delta-V was performed to return the satellite to its targeted geographic 
node. Since longitudinal acceleration always has a constant direction with respect to a specific geographic node, it 
was recognized that targeting for the far linit boundary would Increase the tine between station nalntenance naneuvers. 
Figure X pictures the difference between these two targeting scheaes. But there was one conplicatlon. Due to the 
configuration of the first generation GPS spacecraft, Delta-V's had to occur within specific beta-angle windows. This 
led into research directed at characterizing the perfornance of the i. 1 lb thrusters used for station nalntenance 
maneuvers in GPS. Previous Delta-V planning tools nodeled the thrust produced by the §.l lb rocket engine as an Initial 
thrust followed by a linear decrease in force. Data fron the nanufacturer of the rocket engines showed a thrust build 
up to a peak followed by a decreasing force. The difference between these two thrust models is shown in Figure XI. 

Since small errors in the force model of a Delta-V can translate into large position errors several months down the 
road, a program was written incorporating this new thruster model. Aiming for maximum time between station^keeping 
maneuvers with accurate beta-angle window placement reduced the number of Delta-V's and the total "Unhealthy time due 
to station keeping. 
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MONTH OF YEAR 


FIGURE X: PictoraX Representation Showing The Difference Between the Old and New Methods of Station Maintenance 
Manuever Planning. For The New Method To Work, Accurate Knowledge of Thruster Characteristics is 
Needed To Precisely Target a Beta-Angle Window. 


THRUSTER MODEL 


165 PSl DATA 



TIME FROM IGNITION, SECONDS 
■ MODEL + ACTUAL DATA 

LINEAR MODEL 


FIGURE XI: The Linear Thrust Model Compared to a Model Based On The Manufacturer ' s Test Data for The NAVSTAR 
Spacecraft 0.1 lb Hydrazine Thrusters. Incorporat In* The New Model Into The Planning Process 
for Station Maintenance Maneuvers has Improved Performance. 



The JICS defines total Delta-V efficiency as the observed change in orbit semi -major axis compared to the expected 
change in the semi -major axis. The first generation, research and development (Block I), Mavstar spacecraft tend to 
have Delta-V efficiencies of 92%. The second generation, operational (Block II) satellites generally have Delta-V 
efficiencies of 99%. Many factors can affect Delta-V efficiency including rocket exhaust body impingement, variable 
thruster efficiencies due to temperature, catalyst, engine misalignment differences, or attitude notion and attitude 
thruster firings. Figure XII shows a Block II satellite with aajor axis body reference frame. The 9.1 lb thrusters are 
located on a thruster pod as indicated in this figure. Block I satellites have a similar configuration. One major 
difference between the two types of vehicles Is that the Block I spacecraft only have thrusters located along the Y body 
axis. Block II spacecraft have thrusters aligned along both the Y and X body axes. Given that the Y axis rocket engine 
exhaust Impinges against the satellite body, and Block II Delta-V s are performed using the X body axis thrusters, the 
majority of the Block I Delta-V efficiency loss can probably be attributed to the thrust impingement. The Delta-V 
efficiency has been trended for each satellite and the KCS compensates for Inefficiencies during the planning of a 
Delta-V. Unfortunately, Delta-V efficiency doesn't only vary from vehicle to vehicle, but between Delta-V's as well. 

Delta-V maneuvers are modeled In the Kalman Filter’s reference trajectory for the satellite undergoing the 
maneuver. Before 1989, an impulse model was used to account for the Delta-V in the reference trajectory. The impulse 
model instantaneously changes the reference trajectory at the midpoint time of the orbital maneuver. Though the Impulse 
model was adequate for short duration maneuvers, Delta-V's of varying, and greater lengths than a few minutes were not 
well modeled. A thrust model was implemented which integrates the force of the rocket engines over the duration of the 
bum. This thrust model still uses the old linear force model mentioned above, but the prediction accuracy was still 
better than the impulse model. The thrust model had another advantage; the effects of attitude motion could be 
approximated with the available thruster misalignment Inputs. When rocket engine characterization work began, a 
database of attitude notion during maneuvers also was started. This data could now be incorporated Into the thrust 
model to better improve its prediction accuracy. 
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Until very recently, immediately after a thruster firing, a "Q-Bump" on the satellite's ephemeris states was used 
to account for the errors in the Delta-V model. The parameters used in the covariance matrix for the "Q-Bump" were 
excessively large given the available thrust nodel. Work had begun on optimizing the values to set the covariance 
matrix after a Delta-V until it was discovered that the process noise expression in Equation 14 could be modified with 
values based on the Delta-V duration and expected error. Between the new thrust model, approximation of attitude 
errors, and the addition of maneuver process noise, large Delta-V s are beginning to look like the attitude thruster 
firing in Figure X. Before 1989, satellites were unavailable to the user community for up to four days after a 
maneuver. Today the same maneuver will cause under six hours of unavailability. 

But improvements can still be made. An effort is under way to include the same rocket engine thrust curve pictured 
in Figure XI and used with maneuver planning in the reference trajectory thrust model. A rather simple relation is used 
currently to calculate the amount of process noise. Though this relation performs well, determination whether it's 
optimal needs to be made. In the near future, attempts to approximate the attitude thruster firing which occurs during 

a typical Delta-V will be made. 


CONCLUSION AND RECOMMENDATIONS 

Overall The MCS Kalman Filter is quite robust. The bottom line is that highly accurate navigation data is 
routinely transmitted to a worldwide community of military and civilian users with minimal interruption In service. 

The online monitoring and control of the Kalman Filter is relatively simple. However, some improvements in the 
operation of the MCS^alman Filter could be made. For example, many parameters which reside in the MCS database such 
as the [P] matrix value for "Q-Bumps" or Delta-V process noise would yield more control if they were updatable on the 
online system. The aging of the reference trajectory would be less and ephemerls upload prediction would be better if 
more orbital perturbations were modeled. Better clock predictions should be Possible. 

The MCS should have a capability to monitor user accuracy that is independent of the Kalman Filter. The way ERD s 
and ORD's are calculated assumes that the Kalman Filter solution is truthful. It's easily possible for the MCS 
operators to see large ERD's and ORD's on a satellite which is a result of Kalman Filter corruption. Thanks to the MCS 
partitioning scheme this error can be caught rather quickly as the corruption starts to affect other states in the 
partition. Unfortunately it’s still possible to be fooled and mistakenly upload satellites with corrupted data, 
therefore an independent check is needed. 

For Delta-V' s, it's apparent with the proper engineering data on a vehicle system and associated rocket engines, 
highly accurate trajectory predictions are possible. Though relatively late for the GPS system, satellite attitude 
motion and attitude thruster firings during a maneuver could have been accounted for real time with data from rate 
gyroscopes on board the satellite. This data could have been used in a specialized Kalman Filter routine designed for 
maneuvering spacecraft. To counter the variability in Delta-V efficiencies, it is highly reco-ended to have 
thrusters located on a satellite such that their thrust is directed outwards, minimizing any potential body 

impingement. ^ should be careful not to oversimplify the data available to the operator of systems 
incorporating a Kalman Filter, while not making the operation of the Filter overly complex. The full potential of 
system should also be realized given the desip constraints if at all possible. I believe the MCS Kalman Filter 
experience shows what operations personnel can pull off given the relative freedom to make improvements. 
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